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ABSTRACT 

Cometary ultracompact Hn regions have been modelled as the interaction of the 
hypersonic wind from a moving star with the molecular cloud which surrounds the star. 
We here show that a similar morphology can ensue even if the star is stationary with 
respect to the cloud material. We assume that the H n region is within a stellar wind 
bubble which is strongly mass loaded: the cometary shape results from a gradient in 
the distribution of mass loading sources. This model circumvents problems associated 
with the necessarily high spatial velocities of stars in the moving star models. 
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1 INTRODUCTION 

The ultracompact H n regions (UCHuR) found deep within 
molecular clouds provide important information on the early 
phases of the interaction of massive stars with their natal 
environment. The disruption of the cloud material by the 
hypersonic winds and UV radiation fields of these stars is a 
severe barrier to an understanding of the process of massive 
star formation - only by studying the disruption process will 
anything be learnt about the innermost regions of the pro- 
tostellar cloud. The disruption process also involves many 
important problems of gas dynamics. Considerable theoret- 
ical effort has gone into modelling the varied morphology of 
UCHuR. 

Most theoretical attention has so far been given to the 
cometary regions, which comprise about 20 per cent of ob- 
served UCHuR (Churchwell 1990). Perhaps the most de- 
tailed model, that of Van Buren & Mac Low (1992, and 
references therein), treats them as the steady-state partially 
ionized structures behind bow shocks driven by the winds 
of stars moving through molecular cloud material. Although 
it has been argued that some morphologies which are not 
apparently cometary (in particular, core-halo) can be ex- 
plained as cometary structures viewed close to their axes 
( Mac Low et al. 1991 ), the shell and multiply peaked mor- 
phologies cannot, and so other models should also be inves- 
tigated. 

There are also a number of unresolved questions with 
regard to the cometary models. First, the star is assumed to 
be moving through relatively homogeneous molecular cloud 
gas. Yet it is well known that cloud material has a clumpy 
distribution down to very small scales. Secondly, current 



cometary models require rather high stellar velocities, char- 
acteristically 10-20 km s _1 (Van Buren & Mac Low 1992; 
Garay, Lizano & Gomez 1994). The origin of such a high 
velocity dispersion is proble matic. Moreover, as emphasized 
by Gaume, Fey & Claussen ( 199^ ), if these high stellar ve- 
locities are correct, the vast majority of UCHuR should be 
well isolated from other structures for most of their lifetime. 
Yet they often appear closely associate d with other radi o 
continuum sources and molecular cores (Garay et al. 1994). 



Churchwell ( 1995 ) a rgues that the parameters used by 
Van Buren & Mac Low (1992) have a fairly large range of 
possible values. In particular, the ambient density used in 
these models is significant ly smaller than has be en observed 
ahead of some UCHuR ( |Cesaroni et al. 1994 ). Since the 



shape of the bow shock is determined by the ram pressure 
of the impinging molecular material, such high densities will 
allow cometary regions to form for substantially smaller rel- 
ative velocities. However, such high densities are not com- 
mon, and are likely to be confined to relatively small clumps. 
Since these clumps will not be immediately destroyed on 
their passage through the bow shock, the model discussed 
in the present paper may in fact be a more appropriate de- 
scription of the structure of these regions. 

Finally, as more detailed observations become available, 
the similarity of sha pe to the structures predicted by Van 
Buren & Mac Low (1995) becomes less convincing. Some 



of the cometary UCHuR observed by Gaume et al. (199 



have tails that hook back around the centre of the nebula, 
while the bow shock model predicts structures with bright 
arcs and diffuse tails - Gaume et al. indeed propose that 
such arc-like regions be considered a separate morpholog- 
ical class. When observed at higher resolution, apparently 
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cometary structures often seem to dissolve into collections 
of individual emission knots (e.g. Kurtz, Churchwell & Wood 
1994), rathe r than remaining as smooth st ructu res. 

Dyson (1994) and Lizano & Canto (1995) have sug- 
gested that the interaction of a massive star with clumpy 
molecular material provides a natural alternative model for 
UCHllR, which circumvents the problem of the relatively 
long (10 s yr) lifetimes req uired by their hig h frequency com- 



pared to field OB stars ( Churchwell 1990| ). Clumps act as 
localized reservoirs of gas which can be injected into their 
surroundings by photoionization and/or hydrodynamic ab- 
lation. The ionized region is bounded by a recombination 
front; that is to say, an Hn region is held at a constant ra- 
dius by the high gas density which results from the mass 
loading. The frictional heating of the flow in mass loading 
regions may explain the very high temperatures i nferred for 
some fr acti on of the gas in the Sgr B2 F complex (Mehringer 
et al. lQ. 

There are many possible variant s of this basic premise. 
Dyson, Williams & Redman ( Paper 1 , 1995 - hereafter) have 
described models in which a strong stellar wind mass loads 
from clumps. In the model presented in Paper I (as in the 
present paper), the ionized gas flo w is every where super- 
sonic. Redman, Williams & Dyson (Paper II, 1996 - here- 
after) have discussed models where either the wind is so 
weak or the mass loading is so great that the interior flow 
shocks or is totally subsonic. 

The models of Papers I and II are spherically symmet- 
ric, and reproduce well many observed morphologies (e.g. 
shell-like and centre-brightened). They cannot, of course, 
produce cometary morphologies. We here describe how a 
simple geometrical distribution of mass loading centres can 
generate such morphologies. This model has two important 
properties. First, it again exploits the clumpy nature of 
clouds. Secondly, the star can be at rest with respect to 
the cloud. The structures predicted are dependent only in 
detail on the exact form of the mass loading distribution. 
We consider here only the simplest case, namely where the 
flow in the UCHllR remains supersonic, and defer a discus- 
sion of the much more complex subsonic flows, and those in 
which an initially supersonic wind goes through a termina- 
tion shock, to a later paper. 



2 MODEL 

We adapt the supersonic, dust-free, isothermal flow model 
of Paper I to include a non-spherical distribution of mass 
loading sources. As in Paper I, we take the stellar wind mo- 
mentum output rate, /i*, as 



/J* 



.22 

47rr pv 



const, 



(1) 



where r is the distance from the star, p is the gas density 
and v is the velocity of the flow. In this supersonic model, 
we assume that the streamlines are essentially radial. Hence, 
we assume equation ([[]) holds within any solid angle element 
and take the mass conservation equation to be 



1_ d_ 

r 2 dr 



(r pv) 



qo exp(r cos#/A), 



(2) 



where q is the mass loading rate and A is the scale height 
of the mass loading distribution. In this model, p and v are 



thus functions both of the distance from the central ionizing 
source, r, and of the azimuthal angle, 6. 

The distribution of mass loading has been chosen to 
vary as a function of the distance, x = r cos 8, perpendicu- 
lar to a certain plane through the star (this is not the con- 
ventional correspondence, as the Cartesian co-ordinates are 
oriented by the line of sight assumed to the region, while the 
spherical polars are oriented by its symmetry axis). At large 
positive x, the mass loading is great, whereas it decreases 
to a negligible level at large negative x. Such a distribution 
of mass loading might be expected where a young ionizing 
star formed at the border of a molecular cloud. Closer in to 
the centre of the cloud, an increased density of mass loading 
clumps will produce a higher mean mass loading rate; away 
from the cloud, the mean density of material will decrease 
and the mass loading rate will drop. 

Integrating equation (g), we find that, once the mass 
flux is dominated by the mass loading material, the flow 
density is 

where r = rcos#/A. The recombination front radius, Rn 
(which is a function of the angle, 9), follows from the 
Stromgren relationship in a solid angle element dfl: 
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where f3 is the case B recombination coefficient, n = p/rh is 
the number density of nucleons (we assume the gas inside 
the recombination front is almost completely ionized), and 
S 1 * is the total stellar output of Lyman continuum photons. 
Hence 



, / W/3# V /U Tl ,i/ii 

cose = x [s^J /(tr) 

A 



(•>) 



8xl0 16 cm 



\ .4/11 tf-i/ii ■ -2/11 T , a/n fc \ 

) 9-30 S 48 M28 J ( r R) I I 6 ) 



where tr = _Rr cos 0/ A and 



7(r) 



[(t 2 -2t + 2) exp(r)-2] 



4 dr 



(7) 



In equation (6), we have taken (3 — 2x10 13 cm 3 s 1 , m 
2x 10~ 24 g, and defined scaling variables 5*48 = SV/10 48 s" 



q/W 



30 gcm- 3 s- 1 (cf. 



/t28 = A*/10 gems - and 
Paper I). 

In Fig. we plot J(r) 1 ^ 11 as a function of r. From equa- 
tion (5), the recombination front directly into and against 
the mass loading gradient occurs where J(r) = ± const: this 
will be symmetric if tr <C 1, but, as tr increases, the ion- 
ized region will soon become highly asymmetric. The tail 
of the curve at large negative r tends to a constant value 
of I(t): for high enough ionizing fluxes, the flow in the di- 
rection of low mass loading is essentially free and the gas 
remains ionized at all radii (within the assumptions of this 
model) . 

In Fig. B we plot the loci of Rr(9) (i.e. the envelopes 
of the UCHllR when seen from a direction perpendicular 
to the density gradient), for a range of UCHllR sizes rela- 
tive to the exponential scale height - these curves are given 
parametrically by equation (5). 
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Figure 3. Grey-scale plot of emission from UCHnR with non-uniform mass loading, with the recombination front in the direction of 
least loading (9 = 7r) at (a) Tt a ii = —0.3, (b) r ta ;i = —1. The star is at the intersection of the axes. Contours are also shown, at 2.5, 5, 
7.5, 10, 20, 30, 40, 50, 60, 70, 80 and 90 per cent of the peak intensity [the faintest contour away from the edge of (a) is at 50 per cent, 
in (b) at 30 per cent]. 




Figure 4. Grey-scale and contour emission measure image of a cometary model, for r = —5 in the 'tail'. The contours are chosen as 
in Fig. ^. To the right, line profiles are shown, calculated at the various labelled points. The thin curves are the zero-broadening line 
profiles, while the bold curves are shown for thermal broadening appropriate to a radial Mach number of 2.2 at the recombination front 
at the head of the cometary shape. The velocity scale is constant between the profiles, but the peak intensity scale is arbitrary — the 
total strength of the emission for each observation may be judged from the emission measure value at that point. 
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Figure 1. The function /(t) 1 / 11 , which determines the shape of 
the cometary regions (cf. equation (5)). 




Figure 2. Cometary shapes, for a range of UCHuR sizes relative 
to the scale height of the mass loading distribution. The star is, 
in all cases, at the intersection of the axes. 



2.1 Emission measure 

In this paper, we consider only the simplest case of cometary 
structures observed edge-on. We define Cartesian coordi- 
nates (x, y, z) with the star at the origin, where (x, y) are in 
the plane of the sky and x = rcosf?. The emission measure 
at x = bi, y = 62 is given by the integral along the line of 
sight (parallel to z) 



£(61,62) = 



.4 



n dz 



i-iv 



(t 2 - 2r + 2) exp(r) - 2 



/ r 8 dz, 



(8) 



(9) 



where r = 61 /A and zr is given by z R = b\ tan 2 #r — 6§, with 
#r given by equation (5). 

Example plots of the emission measure are shown in 
Figs ^| and ^, with contours chosen to be similar to those used 
in observational papers. The size of the cometary structures 
is parametrized by run, the value of tr in the direction of the 
tail, 9 — tt. The two cases shown in Fig. [] demonstrate the 
appearance of the cometary form as the size of the UCH11R 
increases relative to the scale of the non-uniformity in mass 
loading. For r ta ii = —0.3, the region is mildly aspherical 



and may be compared to the models of Paper I. For r ta ii = 
— 1, the structure is similar t o tha t of the arc-like UCH11R 
described by Gaume et al. (1995). When r ta ii = —5, the 
morphology is cometary. 



2.2 Line profiles 

Optically thin line profiles can be calculated from the equa- 
tion 

£(&i,& 2 ,«i) = / n 2 f f^v r - Vi\ dz, (10) 

J— Z R 

where (&i , 62) are the coordinates (x, y) of the point observed 
in the nebula, and the line-of-sight component of the velocity 
is v\. The distance r is given, as above, as a function of z. The 
flow velocity, v r (61,62,2), is given by equations ([j]) and ^ 

as 



r 3 4nqo 



6? 



(& 2 -2bi+2)exp(6i)-2 



(11) 



i.e., v r = A(bi)/r s , where A(bi) is a constant for each line 
of sight. 

In the limit of zero broadening (where f(v) = 5(v), 
the Kronecker delta function) , contributions to the integral, 
equation (p"c|), come at solutions, Zi, of the quartic equation 

(b 2 + zlfvi = Zi A{bi), (12) 

in the range [— zr, zr], where b 2 = 6f + 62. The line profile 
is then given by 

(b 2 + z 2 ) 3 



S(bi,b2,vi) = [71(61,62,; 



A(6i)|& 2 -3z 2 



(13) 



Line profiles including broadening can be obtained by 
convolving this line profile by a line-of-sight velocity function 
g(v). This integration will remove the singularity in equa- 
tion ( |l3| ) where z 2 = 6 2 /3. For instance, thermal broadening 
is given by convolving with a Gaussian 



g(v) = 



2iva 



exp( 



2 /2c 



(14) 



where the one-dimensional velocity dispersion, a, is given by 
a 2 = ksT/Mz = fhc 2 /Mz for line emission from a species of 
atomic mass Mz, if fh is the mean mass per particle and c is 
the isothermal sound speed. If the integral, equation (hp ) , is 
performed directly with f(v) = g(v) given by equation (14j), 
major contributions to the total will come from near the 
solutions Zi of equation (^) and from line-of-sight velocity 
extrema. 

In Fig. ^, we present line profiles for a selection of points 
within the cometary shape. Close to the edge of the region 
(points A, H), the flow is nearly in the plane of the sky, and 
the line widths are only slightly above the thermal widths. 
For lines of sight further into the nebula (points B, C, D), a 
larger component of the flow velocity will be seen: the line 
profiles become increasingly double-peaked. Further around 
the bright shell (point E), the double-peaked profiles are sim- 
ilar to those found along a line to the apex. The unbroadened 
profiles at points C and F show the sharp peak at the max- 
imum line-of-sight flow velocity, but when smoothed by the 
thermal broadening this has little obvious effect. At points D 
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Figure 5. Channel maps for cometary UCHiiR, shown for velocities v as fractions of the radial velocity at the recombination front, 
v r jj, at the head of the comet (8 = 0). The profiles are thermally broadened assuming that this velocity is 22 km s —1 , i.e. Mach 2.2. The 
grey-scale representation is normalized in each image; the contours are scaled to 2.5, 5, 7.5, 10, 20, 30, 40, 50, 60, 70, 80 and 90 per 
cent of the peak intensity in the v = channel map. The ticks on the axis are at intervals of A, the scale height of the mass loading 
distribution; the dashed contour at the edge of the region is as in Fig. M. 



and G, this spike is outside the range of velocities displayed 
in Fig. ^: the line profile has two fairly well-separated com- 
ponents, with noticeably stronger tails to high velocity. 

In Fig. we present the line profile data in the form of 
maps at a single velocity (similar to observational 'channel 
maps'). These are shown for single velocities, v, but, since 
the emission has been thermally broadened, they show the 
structures expected for narrow-bandpass observations. At 
low velocities, the emission is bright, and is dominated by 
the dense gas at the edge of the UCHiiR. As the relative 
velocity of the bandpass increases, however, the emission 
peak gradually draws away from the edge of the region and 
fills in the central hollow. 

This displacement of the peak of the diffuse emission 
with increasing relative velocity is a characteristic of radial 
flow models such as that discussed here. In contrast, in the 
bow shock and 'champagne' models, significant line-of-sight 



velocities can be produced at the edge of the region by flows 
around the shell structure. The finite spatial resolution of 
radio observations and the presence of strong, discrete emis- 
sion components will, however, confuse this as an observa- 
tional diagnostic. As can be seen from Fig. ^, the emission 
will only begin to become centre-brightened at relative ve- 
locities <; 22kms~ 1 , and with emission intensities 10 per 
cent of the peak in the zero- velocity image. 

These conclusions are, however, strongly dependent on 
the model assumption that the dominant mechanism which 
broadens the emission is thermal broadening. As we argued 
in Paper I, the emission from supersonic mass loaded winds 
is likely to be dominated by the dense gas near to individual 
clumps. The effect of the finite number of individual clumps 
will be more noticeable in velocity-resolved observations of 
UCHiiR. Also, the tails are likely to be very dynamically 
active, with characteristic turbulent velocities similar to the 
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shear velocity between the wind flow and the (small) clump 
velocity. A dependence of this type would serve to fill in, 
at least in part, the faint central regions of the low- velocity 
maps. Firm conclusions must, therefore, await a detailed 
treatment of the line emission from the regions around in- 
dividual mass loading clumps, which is beyond the scope of 
this paper. 



3 CONCLUSIONS 

In this paper, we have shown how cometary morpholo- 
gies may readily be produced in a mass loaded model for 
UCHllR, with no need for relative motion between the cen- 
tral young OB star and the molecular cloud material, if 
mass loading sources are concentrated to one side. Veloc- 
ity channel maps show a shell structure stronger than the 
integrated emission at small velocities relative to systemic, 
becoming centre-brightened at higher relative velocities. We 
have assumed that the lines are broadened only to the ther- 
mal width: in reality, however, the velocity structure may be 
dominated by turbulence in the dense gas around individual 
mass loading sources. 

We have assumed an exponential distribution in the rate 
of mass loading with distance, x, away from a plane contain- 
ing the star. The general shape of the regions is not strongly 
dependent on the detailed form of the mass loading distri- 
bution, so long as it has a general decrease in one direction. 
In one direction, the recombination front is trapped close to 
the star by a high rate of mass loading, while in the oppos- 
ing direction the material remains ionized to a rather larger 
radius, or even out to infinity. Such distributions might re- 
sult not only near the edge of a molecular core, but also if a 
small relative motion carried clumps into the UCHllR which 
were then evaporated as they crossed it. 

Related models can be considered for different distribu- 
tions of mass loading centres: for instance, if the luminous 
star formed on a ridge of clumpy molecular material (or sur- 
rounded by a thick disc) , the resultant UCHllR would have 
a bipolar morphology. This can be pictured by imagining a 
mirror placed at the plane x = of Fig. ^, to produce a 
structure of two back-to-back cometary tails. 

The model discussed here will give rise to observable ve- 
locity gradients in the head-to-tail direction when the sym- 
metry axis is inclined to the plane of the sky. Since the ve- 
locity of gas crossing the recombination front, wr, is nearly 
constant around the comet (to within 50 per cent even when 
the tail is 50 scale lengths away from the star), the maxi- 
mum observed difference in mean line-of-sight velocit y wil l 
be ~ 2«Rsini, for inclination angle i. Garay et al. (1994) 
find differences of 8-12 km s 



between the mean velocities 
of emission in the head and tail of three cometary UCHllR, 
which can be fitted with a mild inclination of the cometary 
axis to the plane of the sky. 



While Garay et al. (1994) find no velocity gradients 



across the regions they observe, Gaume et al. (1994) list 



three well-defined cometary regions which have an asymmet- 
ric velocity gradient perpendicular to the continuum sym- 
metry axis of the UCHllR. This observation is problematic 
for the model discussed here, as it is also for the bow shock 
model. These contradictory results may relate to morphol- 
ogy: Gaume et al.'s sample tend to have significantly longer 



and more complex tails, while Garay et al.'s sample are more 
arc-like. 

We find that such arc-like morphologies are produced 
for milder asymmetries in mass loading. They form a natu- 
ral intermediate between shell and cometary morphologies. 
Indeed, we would expect arc-like morphologies also to be 
found in bow shock models for UCHllR, if the relative ve- 
locity of cloud and star were low: such a model might, how- 
ever, have too low a gas pressure to correspond to observed 
UCHllR. 



Kurtz et al. ( 1994 ) find some observational evidence 
for a size-density relationship for spherical and unresolved 
regions, but not for cometary and core-halo UCHllR. This 
suggests that a Stromgren relation might hold for the former 
class, but not for the latter. However, given the sparse data 
on which these correlations are based (which do not allow, 
for instance, for correction for the luminosity function of the 
stars) , we do not view this as a conclusive argument against 
the model discussed in the current paper. 

In subsequent papers, we will address the effects of tan- 
gential pressure forces in mildly supersonic models and the 
structures expected when the anisotropic mass loading is 
sufficiently strong to produce a subsonic region of the flow. 
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